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ABSTRACT 

In this paper, we check in detail the validity of the widely adopted point source approxima- 
tion for nearby cosmic-ray (CR) sources. Under an energy independent escape model for CRs 
from the sources, we show that for young nearby sources, the point source approximation 
breaks down at lower energies and the CR spectrum depends on the size and the morphology 
of the source. When applied to the nearby supernova remnants (SNRs), we find that the ap- 
proximation breaks down for some of the individual remnants like the Vela, but interestingly 
it still holds good for their combined total spectrum at the Earth. Moreover, we also find that 
the results obtained under this simple approximation are quite different from those calculated 
under an energy dependent escape model which is favored by diffusive shock acceleration 
models inside SNRs. Our study suggests that if SNRs are the main sources of CRs in our 
Galaxy, then the commonly adopted point source model (with an energy independent escape 
scenario) appears flawed for CR studies from the nearby SNRs. 
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1 INTRODUCTION 

Cosmic rays (CRs) with energies below the knee region (~ 3 x 
10^^ eV) are considered to be of galactic origin. Although the exact 
nature of their sources are not known, the most favorable candidates 
are the supernova remnants (SNRs). They are known to occur in our 
Galaxy at the rate of ~ 1/30 yr^^ with each explosion releasing 
a total kinetic energy of ~ 10^^ ergs. If approximately 10% of 
this energy is converted into CRs, then the total power release is 
sufficient to maintain the CR energy density in our Galaxy which 
is measured to be around 1 eV cm^"^. 

It is also now theoretically established that SNR shock waves 
can accelerate CRs up to energies close to the knee by diffusive 
shock acceleration (DSA) mechanism (Bell 1978, Blandford & 
Eichler 1987). In a simple planar shocks model, such a mecha- 
nism naturally leads to a power-law spectrum of the form E^^ 
with the exponent F = 2 for strong shocks. This value is found to 
be in good agreement with the radio observations of SNRs (Green 
2009). In addition, direct evidence for the presence of high energy 
particles up to few TeVs (1 TeV= 10^^ eV) inside SNRs comes 
from the detections of non-thermal X-rays (see e.g., Parizot et al. 
2006, Bamba et al. 2006) and high energy TeV 7-rays from several 
SNRs (e.g., Aharonian et al. 2006, 2008a, 2008b, 2008d, Albert 
et al. 2007). The non-thermal X-rays are best explained as syn- 
chrotron emission from high energy electrons while the origin of 
the TeV 7-rays is still not certain between the leptonic (via inverse 
compton process) and the hadronic scenarios (via tt'' decays). If the 
high energy 7-rays are of hadronic origin as indicated by the recent 
observations of several SNRs by the FERMI experiment (Abdo et 
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al. 2009, 2010a, b, c), then the measured 7-rays can provide di- 
rect informations about the spectral shape of the primary particles. 
But, TeV measurements made by the new generation Cherenkov 
telescopes like the HESS, MAGIC and VERITAS have found that 
many SNRs show F ~ (2.3 — 2.7) which is steeper than the ex- 
pectations from DSA theory. The discrepancy becomes even more 
severe if we compare with the results of non-linear DSA theory 
which predicts a spectrum flatter than F = 2 (Caprioli et al. 2010 
and references therein). Although this discrepancy is still not yet 
fully understood, for our study we will assume that SNRs are the 
main sources of CRs in our Galaxy. 

Quite often, theoretical studies on the propagation of CRs as- 
sume the sources to be stationary and continuously distributed in 
the Galaxy. This simple assumption seems reasonable for calculat- 
ing the Galactic average CR properties and for understanding the 
diffuse radiations produced by the interaction of CRs in the inter- 
stellar medium (ISM). But, for CR studies in the vicinity of the 
sources where the influence of the source is expected to dominate 
over the background produced by the distant sources, the discrete 
nature of the sources (both in space and time) may become im- 
portant. For instance, in the study of gamma-ray emissions from 
the environment of the sources or from molecular clouds associ- 
ated with them, the emission can be strongly dependent on the age 
and the distance of the source as discussed in Aharonian & Atoyan 
1996, Gabici et al. 2009, Casanova et al. 2010. 

Similarly, for the study of CRs observed at the Earth, the uni- 
form source distribution looks proper only for the distant sources 
but for the nearby sources, a more reasonable treatment would be to 
consider the discrete nature of the sources. For CR electrons at few 
TeV energies, such treatment seems even more important because 
of their fast energy loss rate. Electrons with energies greater than 
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~ 1 TeV cannot travel distances more than ~ 1 kpc in the Galaxy 
through diffusive propagation before they lost all their energies. 
Therefore, high energy electrons from distant and old sources may 
not reach the Earth effectively and it is possible that most of the 
TeV electrons that wc observe arc mostly produced by few young 
nearby sources (Shen 1970, Atoyan et al. 1995, Kobayashi et al. 
2004, Delahaye et al. 2010 etc.). Also for the CR protons and other 
nuclear species which do not suffer significant losses (the typical 
nuclear fragmentation loss time scale in our Galaxy ~ 10^ yr) and 
for which we expect a strong background from the distant sources, 
the discrete treatment of the nearby sources can still be important 
especially at higher energies (Strong & Moskalenko 2001, Bushing 
et al. 2005, Erlykin & Wolfendale 2006 and references therein). 
This is because high energy CRs diffuse relatively faster compared 
to the lower energy ones and hence, they are expected to produce 
stronger fluctuations on their observable properties Uke the spec- 
trum and the anisotropy (Thoudam 2008). Moreover, at these ener- 
gies, the contribution from the recent sources may dominate and the 
effect of discreteness in time may also become important (Taillet et 
al. 2004). 

In most of the studies mentioned above, the discrete sources 
are assumed to be point-like, thereby neglecting their finite size 
and the morphology. At first sight, the point source approximation 
seems reasonable for sources whose size s <C r, the distance from 
the Earth. But, for those whose size is comparable to the distance, 
the point source approximation may break down and it looks more 
appropriate to take their size and morphology into account. Under 
the standard DSA theory, CRs are confined within the SNRs due 
to the strong magnetic turbulence generated by the CRs themselves 
and therefore, it is reasonable to assume that CRs remain confined 
as long as the shocks remain strong enough to act as an efficient 
accelerator. For a typical interstellar matter density of 1 H cm~^, 
the confinement last until the SNR age around 10^ yr (Berezhko 
& Volk 2000). In reality, an energy dependent confinement/escape 
scenario is expected (see e.g. Ptuskin & Zirakashvili 2005, Caprioli 
et al. 2009). Using the Sedov relation between the SNR age and 
the shock radius, if we assume an initial shock velocity of 10^ cm 
s~^, we can roughly estimate that at the age of 10^ yr the remnant 
expands to a size of radius around 100 pc. Such a size is comparable 
to the distance of some of the nearest SNRs like the Geminga and 
the Loopl. The distance to the Geminga is estimated to be ~ 157 
pc (Caraveo et al. 1996) and that to the Loopl as ~ 170 pc (Egger 
& Aschenbach 1995). 

The argimient just mentioned is purely based on the geometri- 
cal consideration, i.e., the source size compared to its distance and 
we have not considered any possible effects due to the propagation 
of CRs in the Galaxy. It is now well accepted that CRs undergo 
diffusive motion due to scattering by the magnetic field irregulari- 
ties and the self excited hydromagnetic waves in the ISM. Measure- 
ments of secondary-to-primary (s/p) ratios like the boron-to-carbon 
indicate that the diffusion is energy dependent with the diffusion 
coefficient increasing with energy (see e.g. the experiments listed 
in Stephens & Streitmatter 1998). If we also take into account such 
an energy dependent diffusion, the validity of the point source ap- 
proximation may become somewhat relaxed for the high energy 
particles, i.e, it may still represent a good approximation even for 
the nearby sources at higher energies. We will discuss this in detail 
later in the paper. 

Recently, Ohira et al. 2011 highlighted the importance of the 
finite source size in the study of gamma-ray emission from SNRs 
interacting with molecular clouds. They claimed that the observa- 
tions of different ganmia-ray spectra from four SNRs (W5 IC, W28, 



W44 and IC 443) by the FERMI experiment could be an effect of 
finite size of the SNRs. In this paper, we will investigate the im- 
portance of the source size for the nearby SNRs considering the 
CR spectrum expected at the Earth. Although SNRs can have com- 
plex morphologies, that too different from each other, for simplicity 
we will consider a spherical geometry for our study. In one part, we 
will consider an energy independent escape of CRs from the SNRs. 
This is discussed in sections 2 and 3. In another part of our study, 
we will investigate the energy dependent escape model under which 
CRs of different energies are assumed to escape at different times. 
This study is given in section 4. Then in section 5, we apply our 
study to the nearby known SNRs and compare the results obtained 
under the different source models. Finally in section 6, we present 
an overview of our results and discuss their implications. 



2 CR PROTON SPECTRUM FROM AN SNR 

In the diffusion model, neglecting losses due to interactions in the 
ISM, the propagation of CR protons in the Galaxy can be described 
by (see e.g. Gaisser 1990 and references therein), 

dNp 
dt 



V • (DVTVp) + Qj,= 



(1) 



where Np(r,E,t) is the differential proton density, E is the ki- 
netic energy, D{E) is the diffusion coefficient and Qp(r, E,t) is 
the source term, i.e. the proton production rate from the SNR. In 
Eq. (1), we also neglect other effects which are relevant mostly be- 
low a few GeVs like the convection due to the Galactic wind and the 
reaccelcration by the interstellar turbulence. We assume the diffu- 
sion coefficient to be spatially constant throughout the Galaxy and 
take D{E) = Do{E/EoY for E > Eo, where Do = 2.9 x 10^^ 
cm^ s"\ So = 3 GeV and 5 = 0.6 (Thoudam 2008). These 
values are different from those given by models based on diffu- 
sive reacceleration. For instance, Trotta et al. 2011 give a value of 
Do ~ 8.3 X 10^* cm^ s"^ and 5 ~ 0.3. 

For sources within a distance of ~ 1 kpc from the Earth which 
are also our main interest here, Thoudam 2007 showed that the 
CR spectrum is not much affected by the presence of the Galac- 
tic boundaries. In fact. Do does depend on the boundaries and is 
proportional to the size of our Galactic halo (see e.g. Trotta et al. 
201 1). For our present study, we neglect such dependencies and 
solve Eq. (1) without imposing any boundary conditions. We then 
obtain the well known Green function Gp{r, r' ,t, t'), i.e. the solu- 
tion for a 5-function source term Qp{r, E, t) = 6(r — r')S{t — t') 
as given below. 
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The general solution of Eq. (1) can be then obtained using, 

/oo rt 
dr dt'Gpir,r',t,t')Qp{r',E,t') 
-OO •' — oc 

The source term in Eq. (3) can be written as, 

Qp{r',E,t') = q{r')q{E)q{t') 



(2) 



(3) 



(4) 



where q{E) is the source spectrum, i.e., q{E)dE is the number of 
protons with energy between E and E + dE produced by the SNR. 
For this part of our study, we assume an energy independent escape 
of CRs from the SNR. We will first consider the burst-like injection 
of particles followed later by the continuous injection case. Later 
on, in section 4 we will discuss the energy dependent escape model. 
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If we assume that the burst-like emission of particles happen at 
time to, we can write the temporal source term as q{t') = 5{t'—to). 
Then, using Eqs. (2) & (4) in Eq. (3), we get, 



q{E) 



8 [7r_D(t - to)] 



3/2 



f OCi 


r-(r'-r)n 
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.4I>(t-to). 



?(r') 

(5) 



The proton intensity can be then calculated using the relation 
Ip{E,t) Ri {c/4n)Np{E,t), where c is the velocity of light. To 
make our calculations simpler, hereafter we take r = 0, i.e., we set 
the origin of the coordinate system at the position of the Earth. 



2.1 Point source approximation 

If we assume the SNR to be a point source located at a distance Vs 
from the Earth, we can write 



g(r') = 5{r' - r,) 
Then, the proton density is obtained from Eq. (5) as. 



Np{E,t) 



8lwD{t-to)f/^ 



exp 



4:D{t - to) 



(6) 



(7) 



Eq. (7) represents the most commonly adopted solution for CR 
spectrum from a nearby single source. For high energy particles 
for which the diffusion radius defined as rai ff = ^/4:D{t — to) is 
much larger than the distance to the SNR rs, the exponential term 
in Eq. (7) tends to 1 which implies. 
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(8) 



For a power-law source spectrum given by q-p{E) = kpE ^ and 
for D[E) oc E^ , Eq. (8) shows that the spectrum of high energy 
protons reaching us follows Np{E) oc i5~(^+t''). Particles with 
rdiff > rs are those which have already passed the Earth. Those 
with rdiff < Ts are the ones which have not yet reached effec- 
tively due to their slower diffusion and therefore, their intensity is 
comparatively much suppressed. 



2.2 Spherical solid source 

Most of the SNRs are observed to roughly follow a spherical geom- 
etry and they come under three main categories: shell-type, plerion- 
type and composite-type. Shell-type SNRs show bright shell struc- 
ture which expands into the ISM with velocities of ~ (3— 10) x 10* 
cm s~^ (e.g. Cassiopeia A). Plerions also known as pulsar wind 
nebulae have filled center normally a pulsar powering high en- 
ergy particles into the ISM (e.g. Crab Nebula). They do not show 
any shell-like features. Composite SNRs have both shell structure 
and filled center (e.g. IC443). The surface brightness of shell-type 
SNRs in radio as well as in X-rays are observed to peak near the 
surface whereas in plerions, it tends to increase towards the cen- 
ter. We can expect that the high energy particles responsible for 
the radio and the X-ray emissions also follow a similar distribution 
within the renmant. 

Let us now consider a spherical solid source. We believe that 
this source model roughly represents the plerions and the composite 
type SNRs. For this model, if r^, denotes the position of the center 
of the SNR from the Earth and ro represents the position of the 
source CRs with respect to the SNR center, we can write r' in Eq. 



(5) as r' = rs + ro and then rewrite Eq. (5) as. 
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In Eq. (9), the integral over the volume element in spherical geom- 
etry is given by. 
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where R denotes the radius of the SNR. We take the source spec- 
trum in this case as q{E) — qp{E)/V where qp{E) is the source 
spectnmi we took in the case of the point source approximation 
(section 2.1) and V = ^ttR^ represents the total SNR volume. We 
assume that CRs are uniformly distributed throughout the SNR vol- 
ume before releasing into the ISM and we take the spatial source 
term as, 



q{ro) 



1, for ro ^ R 
0, otherwise 



(10) 



Integrating Eq. (9) over and 00, we get 

qp{E) 
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Using the properties sinh(a;) ~ x and — > 1 for very small x, it 
is easy to check that for very small R, Eq. (11) reduces to the point 
source solution (Eq. 7) at all energies. 



2.3 Splierical surface source 

If we assume that the CRs are distributed uniformly only on the 
surface of the SNR before they are released, the spatial source term 
in Eq. (9) can be written as. 



g(ro) = i5(ro - R) 



(12) 



and the source spectrum as q{E) — qp{E)/A, where A — 4itR^ 
denotes the total surface area of the SNR. The CR density in this 
case is then obtained as. 
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zexp 
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(13) 



Here again, we can notice that Eq. (13) tends towards the point 
source solution for very small values of R. The spherical sur- 
face source model considered here closely represents the shell-type 
SNRs and it is probably more relevant than the solid source model 
for CR studies in our Galaxy. It is because according to the recent 
catalogue of Galactic SNRs, 78% of the total known SNRs are of 
shell-type while the remaining 12% and 4% are of composite and 
plerion types respectively (Green 2009). 

In Fig. 1 top panel, we compare the spectra obtained under the 
different source models for t = 10^ yr and 10'' yr. The calcula- 
tion assumes to = 0, the SNR distance as rs = 0.15 kpc and the 
diffusion constant as Do = 2.9 x lO'^** cm^ s~^. For our present 
illustration, we take the source spectral index as F = 2 which is 
the value predicted by DSA theories inside SNRs. Later on, in sec- 
tion 5 when we apply our study to the nearby known SNRs, we 
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E (GeV) 




E (QeV) 

Figure 1. CR proton spectra from an SNR with distance rs = 0.15 kpc at 
different times t = 10^ yr and 10* yr for different source models: point 
source (solid Une), surface source (dashed line), solid source (dotted line). 
We consider an energy independent burst-like injection of CRs at to =0 
and we take r = 2.0, 5 = 0.6 and Bq = 3 GeV. Top: Do = 2.9 x 10^* 
cm^ s-i. Bottom: Do = 2.9 x 10^^ cm^ s'^. 

will use values which are determined based on the observed CR 
data. In Fig. 1 top panel, we can see that for a given value of t, the 
point source solution (solid line) above some energy Ept agrees 
well with the results of the surface source (dashed line) and the 
solid source (dotted line) models, while below Ept the results are 
quite different. Ept is roughly the energy at which rdiff ~ rs. We 
can check that for S > iSpt for which rd,// > r^,, Eqs. (11) & (13) 
tend towards Eq. (8) which is the asymptotic solution of the point 
source approximation at high energies. The bottom panel shows the 
results for larger value of diffusion constant Do = 2.9 x 10^^ cm^ 
s^^. The only difference between the two sets of results is that Ept 
is shifted towards lower values as Do increases. This shows that 
the point source becomes valid over a broader energy range as Do 
takes larger values. In order to understand the effect of the source 
distance, we show in Fig. 2 the results obtain for = 0.3 kpc 
by keeping all other parameters same as in Fig. 1. On comparing 
the results in Fig. 2 to those in Fig. 1, we can see that apart from 
the scaling down of the flux and the right shifting of Ept due to 
the increased source distance, the differences between the different 
source models also become smaller. This is simply the geometrical 
effect mentioned in section 1, i.e. as the source distance increases, 
the point source approximation becomes more valid. 

These results can be understood as follows. The diffusion ra- 
dius Tdiff = \/ D(t — to), which is the effective distance from 
the SNR travelled by CRs due to diffusive propagation, is a strong 
function of Do and t. The larger the values of Do and/or t, the 
larger is rdiff and the energy Ept which satisfy the condition 




E (GeV) 

Figure 2. Same as in Fig. 1 but for rs = 0.3 kpc. 




Age t (yr) 



Figure 3. Source age t versus distance rs plot for Ept = 3 GeV and 10 
GeV energies. The area below each line represents the (rs , t) parameters 

space where the point source represents a good approximation for energies 
E > Ept. We assume Do = 2.9 X 10^* cm ^ s-\ <5 = 0.6, Dq = 3 GeV 
and t() = for our calculation. The black dots represent the nearby known 
SNRs with distances ^ 500 pc (see Table 1). 



rdiff = Vs decreases. Similarly, for larger source distances Vs, we 
can understand that Ept shifts towards higher values. These results 
show that the point source can remain a valid approximation even 
for the nearby sources as long as the particles satisfy the condition 
rdiff » (R, rs). For rs = 0.15 kpc. Do = 2.9 x 10^* cm^ s'^ 
and to = (Fig.l top panel), we obtain Ept = 1.2 x 10^ GeV and 
57 GeV for t = 10^ yr and 10* yr respectively. The corresponding 
values for rs = 0.3 kpc (Fig. 2 top panel) are found to be 1.2 x 10*^ 
GeV and 575 GeV respectively. For a given age or distance, the 
closer or the older the source is, the lower is the Ept. This is shown 
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in Fig. 3 for Ept = 3 GeV (solid line) and 10 GeV (dashed line). 
The area below each line represents the parameters space in (rs, t) 
where the point source approximation works as a good approxima- 
tion for all energies above the given Ept. In the same figure, the 
black dots represent the nearby known SNRs with distances ^ 500 
pc (see Table 1). We can see that only the Loopl and the Geminga 
remnants Ue below both the lines while the Monogem Ue just above 
the 10 GeV line. The other two SNRs, Vela and G299.9-2.9 are lo- 
cated well above the lines. It is worth mentioning that CR spectrum 
below ~ (1 — 10) GeV are very much likely to be modified by the 
solar modulation and hence, only those above this energy region are 
reliable for estimates of Galactic CR properties. Therefore, as far 
as the CRs of our interests are concerned, the point source approxi- 
mation looks valid only for the Loopl, Geminga and the Monogem 
among the nearest SNRs, while for the others it looks important to 
take their sizes into account in the calculations. It should be noted 
that for distant SNRs like SN185 for which the distances rs » -R, 
the point source will always remain a good approximation indepen- 
dent of their ages. We will show this in detail in section 5. 

The results shown in Figs. 1 and 2 are obtained for the burst- 
Uke injection of particles. Let us now investigate the case of con- 
tinuous injection of particles. For continuous injection for a finite 
time interval from to T with the injection rate qdE), the solution 
of Eq. (1) can be obtained by integrating Eq. (5) over the injection 
time as. 



Np{E,t) = 
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exp 



ADit - 1') 

X g(r' 

(14) 



where tf = mm[t,T] and qc{E) = q{E)/T, q{E) represent- 
ing the source spectrum in the burst-like injection case which for a 
point source is given by q{E) = qp{E). For a point source posi- 
tioned at Fs, Eq. (14) becomes, 

1p{E) 



Np{E,t): 



AivrsDT 



{erf(V5i)-erf(Vir)} (15) 



where. 



X2 



and 



Xl 



ADt 



4D{t - tf) 

For t < T, tf = t which implies X2 = oo. Then, using the 
property of the error function, erf (y^) = 1 for X2 = oo, Eq. (15) 
in this case becomes. 



Np{E,t) 



1p{E) 



{1-erf (ViT)} 



(16) 
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For high-energy particles for which the diffusion radius rdi ff S> 
rs, Xl — > 0, and because erf(y^) — > for a; i 0, the spec- 
trvun given by Eq. (16) follows a power law of the form Np{E) oc 
j^-(T+s) ^j^jgj^ jg flatter than the spectrum we obtain in the burst- 
like injection model. A detailed discussion on this topic can also be 
found in Aharonian & Atoyan 1996 in the study of CR spectrum in 
the vicinity of the sources. 

Vort > T,tf = T and X2 = rl/AD{t - T). For particles 
with large rdiff for which X\ — > 0, we can safely write a;2 ^ 1 as 
{t — T) < t. Then, using the property eri{y/x2) 



2yfx2fn for 

^/x2 <C 1, the particle spectnmi (Eq. 15) in this case reduces to 
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The spectral shape of Eq. (17) follows Np{E) oc E~(^+i^^ which 
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Figure 4. CR proton spectra for different source models from an SNR at 
r = 0.2 kpc at different times t = (10^, 10*, lO'') yr. Top: Burst-like 
injection. Bottom: Continuous injection with T = 10^ yr. Other model 
parameters are the same as in Fig. 1 top panel. 



is the same as in the case of the burst-like injection of particles 
discussed earlier (Eq. 8). 

Similarly, using Eq. (14) we also obtain our results for the 
case of the solid and the surface source models by taking into ac- 
count their proper source terms given by Eqs. (10) and (12) re- 
spectively. The results are plotted in Fig. 4 (bottom panel) for 
t = (10^, 10*, 10®) yr along with the results obtained under the 
burst-like injection model (top panel) for comparison. The calcula- 
tions in Fig. 4 assume rs = 0.2 kpc. Do = 2.9 x 10^* cm^s~^, 
r = 2.0 and T = 10" yr. We can see that at all t's, the effect of 
assuming different source models are similar in both the types of 
injection. As discussed above, we can also see that for <; < T the 
spectra in the case of continuous injection are flatter than those in 
the burst-like injection case while at t ^ T, they exactly follow 
the same behavior as shown by the results at 10® yr. 

A short conclusion that we can draw at this stage of our study 
is that for very old sources (t » T), the effect of choosing differ- 
ent source geometry or different particle injection model is negli- 
gible on the CR spectrum. Therefore, the widely adopted burst-like 
point source approximation remains a good approximation for very 
old nearby sources at all the energies. However, for young nearby 
sources, the spectrum at high energies strongly depends on type of 
the particle injection model and at lower energies, it starts depend- 
ing on the physical size and the geometry of the source irrespective 
of the type of the injection model unless the source is really closed 
to the Earth i.e., only a few pc away as shown in Fig.3. 
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3 HIGH ENERGY ELECTRON SPECTRUM FROM AN 
SNR 

The diffusive propagation of higli energy electrons in the Galaxy 
can be described by the following transport equation, 

8 BN 
V-{DVNe) + -^{b{E)Ne} + Qe = -g^ (18) 

where Ng{E, t) is the density of electrons with kinetic energy E, 
b{E) = — di5/diis theenergy lossrateandQe(i', denotes the 
electron injection rate into the ISM. The Green fiinction of Eq. (18) 
can be obtained as given below (see e.g., Ginzburg & Syrovatskii 
1964, Gratton 1972), 
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8[7rf{E,E')f^ h{E) 



exp 



-(r'-r)2 



4.f{E,E') 

X 5[t' -t + g{E,E')\ 
(19) 



where, 
f{E,E')= f 

J E 



E' 



D{u) 



du 



/■^' 1 

and g(E,E') = / -r-^du 
^ ^ Je Kw) 



IE K") 

For our present study, we assume that the energy loss of the elec- 
trons are due to synchrotron and inverse compton interactions 
which are true mostly for energies E>1Q GeV. We take, 

h{E) = aE"^ (20) 

where a = 1.01 x 10^^^(wph + wb) GeV s^^ and, -Wp]-, and 
wb represent the energy densities in eV cm~^ for the background 
photons and the magnetic field respectively. Eq. (20) assumes that 
the inverse compton scattering of the background photons occur in 
the Thompson regime. 

The general solution of Eq. (18) is given by, 

/DO roc pt 

dr dE' I dt'Ge{r,r',E,E',t,t') 
-oo J E J — oo 

X Qe{r',E,t') 
(21) 

For an energy independent burst-like injection of electrons at time 
to, we take the source term asQe(r',-B',f') — q{r')q{E')S{t' —to) 
where q{E') oc E'^^ denotes the source spectrum. Now, setting 
r = as we did for the protons in section 2 and performing the 
integrals over E' and t', Eq. (21) becomes, 

r— 2 — oo 



Ne{E,t) = 



q{E) 

8(7rC)3/2 



_ E. 
~ El 

X exp 



4C 



dr' 
9(r') 



(22) 



where Et = 1 / (a (t — to)) is the energy at which the energy loss 
time is equal to (t — to). 



C = 



D{E) 
a(l - 5)E 



(23) 



and 5 is the index of the diffusion coefficient. Eq. (22) is valid for 
electrons with energies E < Et. For E > Et, Ne = 0. 



3.1 Point source approximation 

For a point source described by Eq. (6) located at a distance r^, the 
electron spectrum at time t can be obtained from Eq. (22) as given 




E (GeV) 

Figure 5. CR electron spectra for different source models from an SNR 
at rs = 0.2 kpc. We assume wmbr = 0-25 eV cm~^, Wop = 0.5 
eV cm~^ and the ISM magnetic field as 6juG. All other model parameters 
are the same as in Fig. 4. Top: Burst-like injection. Bottom: Continuous 
injection. 



below. 



Ne{E,t) = 



8(7I-C)3/2 



E_ 
El 



exp 



AC 



(24) 



where (E) — k^E^^ is the source spectrum for the point source. 
From Eq. (23), we can see that in the energy region E <^ Et where 
the effect of the energy loss is less important, C — > /)(t — to) and 
Eq. (24) tends towards the point source solution for CR protons (Eq. 
7). Therefore, high energy electrons whose Vdiff S> Vs and E <IC 
Et have spectrum which follows Ne{E) oc E~^^~^2^) which is 
similar to the asymptotic solution of high energy protons (Eq. 8). 



3.2 Spherical solid source 

Following exactly the same procedure as for the protons described 
in the previous section, we obtain the electron spectrum for the 
spherical solid source as. 



Ne{E,t) = 



qe{E) 



X 



Jo 



ro exp 



1 - 



E_ 
El 

4C 



exp 



4C 



sinh 



2C 



drn 



(25) 
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3.3 Spherical surface source 

We also obtain the solution for the spherical surface source as given 
below. 



Ne{E,t) 



qe{E)R 



E_ 
Et 



exp 



AC 



X sinh 



2C 



(26) 



It is easy to check that for very small values of R, the solutions for 
the solid source (Eq. 25) and the surface source (Eq. 26) models 
reduce to the point source solution (Eq. 24). 

The solutions we have obtained above are based on burst-like 
injection of electrons from the SNR. For the case of continuous 
injection, the solutions are given by. 



r — oo 

X / dr'exp 

J oo 



1-^ 

Ei 



AC 



g(r') (27) 



where tf and qc{E), they bear the same definitions as defined in 
the case of protons and E[ = l/{a{t — t')). For a point source at 

Tg, Eq. (27) becomes, 



N, 



{E,t)= [ 
Jo 



8(7rC)3/2T V Ei 



1-4 



X exp 



AC 



(28) 



where qe{E) is the source spectrum we assumed in the case of 
burst-like injection (section 3.1). In Eqs. (27) and (28), C is given 
by Eq. (23) but with Et replaced by E't in this case. Here again, 
we can check that for energies E <^ E't, C ^ D{t — t') and Eq. 
(28) reduces to a solution similar to that of the CR protons (Eq. 
15). Therefore, the same discussions we presented in the previ- 
ous section for the protons under the continuous injection model 
also apply to the electrons. At time t < T, the spectrum of 
high energy electrons with E <^ E't and whose diffusion radii 
rdiff > Vs follow Ne{E) oc £-(^+4) at t > T, they fol- 
low Ne{E) oc £"('^+2'^) which is similar to the result obtained 
in the burst-like injection model (section 3.1). More discussions on 
the different types of electron spectra generated by a CR source un- 
der different particle injection models can also be found in Atoyan 
et al. 1995. 

Using Eq. (27), we also calculate the spectra for the solid and 
the surface sources under the continuous injection model. The re- 
sults are plotted in Fig. 5 (bottom panel) for t — (10^, 10'', 10®) yr 
The top panel shows the results for the case of burst-like injection. 
The calculations assume T = 2.0, 5 = 0.6, Da = 2.9 x 10^* 
cm^ s"\ rs = 0.2 kpc, T = 10^ yr and to = 0. The mag- 
netic field in the ISM is taken as 6/xG (Beck 2001). The total 
energy density of the background photon field is assumed to be 
yjph = wmbr + Wop, where wmbr ~ 0.25 eV cm"^ is the 
energy density of the microwave background and Wop = 0.6 eV 
cm~3 that of the ultraviolet-NIR-optical radiation field. The lat- 
ter is taken from the estimates given in Shibata et al. 201 1 for the 
galactocentric distance of 8.5 kpc which they obtain using the data 
provided by GALPROP (Porter et al. 2008). On comparing Fig. 5 
top and bottom panels, we can notice that stt < T apart from 
the difference in the slope of the spectra, there are sharp spectral 



breaks present in the case of burst-like injection. These are due to 
the effect of fast energy loss rate for high energy electrons. Elec- 
trons with energy E > Et we, lost before reaching the Earth. We 
also notice that the differences between the different source models 
below Ept are similar in both the types of injection model. At very 
late times t^T, the spectra becomes independent of the injection 
model or of the somce model and they exhibit the same shapes with 
breaks atE — Et. 

On comparing the results of electrons (Fig. 5) to those of the 
protons (Fig. 4), we can see that except for the presence of spectral 
breaks in the case of electrons, the results are quite similar in all 
other respects. Even the differences between the results obtained 
for different source models are similar. Therefore, the overall con- 
clusions on the validity of the point source approximation that we 
had drawn earlier for the protons also apply to the electrons. 



4 ENERGY DEPENDENT CR ESCAPE FROM SNRS 

So far, we have only considered a simple model of CR escape 
from the SNRs where CRs of all energies are assimied to escape at 
the same time independent of energy. However, detailed theoretical 
studies suggest that their escape mechanism can be more complex 
which may be strongly related with the acceleration process itself 
and depends on the shock dynamics as well as on the particle en- 
ergies and their back reaction on the shocks (see e.g. Malkov & 
Drury 2001, Ptuskin & Zirakashvili 2005). 

As already mentioned in section 1, under DSA theory CRs 
are assumed to be confined by the magnetic turbulence generated 
by the CRs themselves. They cannot escape the remnant as long 
as their upstream diffusion length normally defined as Idiff = 
Ds{E)/us is less than the escape length from the shock front 
which is usually taken as lesc = (i.Rs, where -u,, and R,, denote 
the shock velocity and the shock radius respectively and the con- 
stant ^ ~ (0.01 - 0.1) (Ptuskin & Zirakashvili 2005 and refer- 
ences therein). In the Bohm diffusion limit where the maximum 
confinement is acheivable, the upstream diffusion coefficient de- 
pends on the particle energy and the upstream magnetic field 
as Dg{E) oc E/Bs. Under this condition, the escape energy E^sc 
follows. 



Eesc OC BsRsUs 



(29) 



There are strong theoretical arguments which suggest that CRs 
might amplify the magnetic fields near the shock surface (see e.g. 
Caprioli et al. 2009). This idea is also supported experimentally 
by the recent observations of thin X-ray filaments inside several 
SNRs, which are most likely synchrotron emissions of high energy 
electrons in the presence of strong magnetic fields of the order of 
~ (100-1000)^0 (Volket al. 2005). Taking such possible ampli- 
fication into account, we can assume that the magnetic field scales 
with the shock velocity as Ba oc uf , with the index d representing 
the degree of amplification. Some studies suggest that d can reach 
values as high as 1.5 (Bell 2004). 

One reasonable assumption of DSA theory is that CRs do not 
escape during the free expansion phase of the SNR evolution. It 
is because shock waves traveling at some constant velocity can al- 
ways overtake particles undergoing diffusive motion (Drury 201 1). 
However, during the Sedov phase when the shock velocity de- 
creases with the age t as oc t^*'" and the shock radius increases 
as Ra oc t^'*', some of the high energy CRs can start escaping 
because of their relatively larger diffusion length {Idiff > lesc)- 
Therefore, imder the Sedov scaling, the escape energy at any stage 
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during the evolution can be obtained using Eq. (29) as. 



This gives, 



Eesc OC t 



-(0.2+0.6d) 



for d = 1.5 



(30) 



(31) 



In deriving Eq. (31), we assume that Ds{E) scales linearly 
with E. But, the exact dependence is still not well understood and 
depends on some poorly known yet important quantities like the 
spectral distribution of the self-excited turbulence waves, their dis- 
sipation rate and their CR scattering efficiencies. Moreover, mag- 
netic field amplification and the dynamical reaction of the acceler- 
ated particles on the shock structure are also not fully understood. 
Due to these uncertainties, a simple but reasonable approach which 
is commonly followed is to parameterize the escape energy as given 
below (Gabici et al. 2009, Ohira et al. 2011), 



t 

^sed 



(32) 



where Emax is the maximum CR energy and tsed denotes the start 
of the Sedov phase. We assume Emax = 10® GeV (= 1 PeV) and 
taed = 500 yr for our study. Eq. (32) assumes that the escape of 
the highest energy particles start at the onset of the sedov phase 
itself. For detailed studies of particle escape from SNRs, see e.g. 
Ptuskin & Zirakashvili 2005, Caprioli et al. 2009, Caprioli et al. 
2010. Using Eq. (32), we can calculate the escape time tesc as a 
function of energy as, 

f E \ 

Uso(E) = tsed -E, (33) 

\ ^max / 

At some later stage of the SNR evolution when the shock 
slows down and does not efficiently accelerate the CRs, the tur- 
bulence level in the vicinity of the shock goes down and no par- 
ticles can remain confined effectively within the remnant. At this 
stage, we can asstmie that all the CRs escape into the ISM. As pre- 
viously mentioned, for an ISM density of 1 H cm^"*, this happens 
at around 10^ yr after the supernova explosion (Berezhko & Volk 
2000). Taking this into account, the CR escape time for our study 
is taken as. 



Tasc{E) = min [tesc(E)., Ufyi] 



(34) 



Using the Sedov relation between the shock radius and the 
SNR age, we can also calculate the escape radius Resc which we 
define as the radius of the SNR at the time when CRs of energy E 
escape as follows. 



Resc{E) = 2.5vo tsed 



Tesc 
tsed 



-0.6 



(35) 



In Eq. (35), uo represents the initial shock velocity, i.e the velocity 
att = tsed which we take as 10^ cm/s for our study. 

Eq. (34) is plotted in Fig. 6 (top panel) where different lines 
correspond to different values of a: solid (0.2), dashed (1.1) and 
dotted (2.0). The plots show that even for the fixed values of Emax 
and tsed, the energy dependence of Tesc{E) strongly depends on 
the value of a. For a = 0.2, except for particles with energies 
greater than 3.5 x 10® GeV all the particles remain confined till the 
end of the SNR evolution. As the value of a increases, lower energy 
particles start escaping at relatively early stages. For a = 1.1 and 
2.0, only particles with energies up to 3 x 10^ GeV and 25 GeV 
respectively are confined till the end of the evolution. The bottom 




E (GeV) 



Figure 6. CR escape time (top) and escape radius (bottom) for an SNR 
under the energy dependent escape model for different values of a: 0.2 
(solid line), 1.1 (dashed Une), 2.0 (dotted line). We assume Emax = 10® 
GeV and tsed = 500 yr. 



panel shows the corresponding values of Resc{E) calculated using 
Eq. (35). For the assumed value of vo, CRs escape starts when the 
remnant expands to a radius of ~ 5 pc and continues until it ex- 
pands up to ~ 100 pc. The latter value denotes the maximum CR 
confinement radius in our study. 

For our calculations in the following, we will assume that at 
the time of escape from the SNRs, CRs are distributed uniformly at 
the shock surface. This assumption is similar to that of the spherical 
surface source discussed in sections 2 and 3. 



4.1 CR proton spectrum from an SNR 

For CR protons, the source term components in the energy depen- 
dent escape model can be written as. 



q{ro) = (5(ro - Resc) 

q{t') = 5{t' - Tesc) 



(36) 



where qp{E) is the point source spectrum given in section 2.1, 
Tesc and Resc are given by Eq. (34) and Eq. (35) respectively and 
Aesc ~ -iirRisc- Now, the proton spectrum in this case is obtained 
using Eq. (13) by substituting the above source parameters as. 



Np{E,t) = 



qp{E) Re 



rsAese\/nD{t-Te^ 



=exp 



esc ~ ' s I 



X sinh 



AD{t - Tesc) 
\2D{t-Tesc), 



(37) 
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Figure 7. CR proton spectra at different times t = (10'', 10'', 10^, 10°) yr 
under the energy dependent escape model for a = 1.1 (soUd Une) and 2.0 
(dashed Une). We assume = 0.2 kpc, r = 2.0, Do = 2.9 x 10^* cm^ 



, 5 = 0.6, Eo = 3 GeV, 



10*= GeV and t^ed = 500 yr. 



In Fig. 7, we show the proton spectra calculated using Eq. 
(37) for a source at = 0.2 kpc at different times t — 
(10^ 10*, 10^ 10'') yr. The dashed lines correspond to a = 1.1 
and the solid lines to a = 2.0. The sharp breaks in the spectra are 
due to the effect of the energy dependent escape time of the par- 
ticles. Particles with energies below the breaks have not yet been 
escaped from the SNR or even if they do, they have not yet reached 
the Earth at the given time t. The effect of choosing different values 
of a is clearly visible. For a = 1.1, the spectra at all times except 
for i 3> 10^ yr peak at relatively higher energies compared to those 
for Q = 2. This is because particles at all energies except for those 
which remain till the end of the evolution are confined for relatively 
longer period in the case of a = 1 . 1 (see Fig. 6 top panel). Looking 
into the individual spectrum, we can also see that at high energies 
it follows a power-law spectrum as Np{E) oc iJ" i . This can 
be understood from Eq. (37) w hich shows t hat for particles with 
large diffusion radius rdi// oc y/ D{t — Tesc), the solution reduces 
to that of the point source approximation at high energies (Eq. 8). 
An additional effect of large rdi ff 'i& that the spectra at the high- 
est energies for the two different a's are very similar. These high 
energy particles are those which escaped the reirmant long ago and 
their Vdifs S> (-Resc, Vs) so that they have already passed the Earth 
at the given time. For these particles, the expected spectrum is al- 
most independent of the chosen values of a, T^sc and Resc- This 
is more clearly visible in the results obtained for t = 10^ yr where 
the two spectra are almost identical to each other over the energy 
range considered here. 

Although taking different values of a result into different 
types of spectrum especially at the lower energies at a given time, 
hereafter we will adopt a = 2.0 for our study. The effects of choos- 
ing other values of a on our results will be discussed later in section 
6. 



4.2 High energy electron spectrum from an SNR 

To proceed, we recall Eq. (26) which represents the electron spec- 
trum Ne{E) for the spherical surface source obtained under the 
energy independent escape model. In that equation, electrons of 
energy E observe at time t had an initial energy E' at the time of 
their escape given by. 



E' = 




E (GeV) 



Figure 8. CR electron spectra at different times under energy dependent 
escape model for a = 2.0. Other model parameters are same as in Fig. 7. 



where to denotes the escape time from the SNR. We can reverse 
the situation and calculate the energy of an electron after time t for 
a given initial energy E' as, 

E' 



E ■■ 



(39) 



l + aE'{t-to) 

For an energy dependent escape, we can now substitute to by 
Tesc{E') and rewrite Eq. (39) as follows, 

E' 



E : 



(40) 



1 + aE' [t - T,sc(E')] 

If qe{E') represents the source spectrum of electrons with initial 
energy E' which escape the remnant at time Tesc{E'), their energy 
E at time t is given by Eq. (40) and their number density is obtained 

using Eq. (26) as given below. 



= [l + aE'{t-TL,)y 



X exp 



AC 



sinh 



r R' 

' s-i-^esc 

2C" 



where C' is given by, 
D{E') 



C = 



-[l-[l + aE'{t-TL.)\'~^} 



(41) 



(42) 



1 - aE{t - to) 



(38) 



a(l - 5)E' 

and flLc = R.sc{E'), = n,,iE') and A',,, = A,,,{E'). 

Using Eq. (41), we calculate the electron spectra at different 
t's for a source distance Fs = 0.2 kpc. The results are shown in Fig. 
8. On comparing with the results obtained for the protons shown in 
Fig. 7 (a = 2.0), one can notice that the major difference is the 
presence of additional breaks at higher energies which are due to 
the effect of radiative energy losses. The breaks at the lower ener- 
gies which are due to the effect of T^sc are seen at the same ener- 
gies for both the type of particles. The electron spectrvun between 
the breaks follow an exponent F -I- |(5 similar to the proton spec- 
trum and at very late times (say at t = 10® yr), it also tends towards 
the point source solution. These results show that also in the case of 
energy dependent escape scenario, for very old sources (t ^ 10^ 
yr) the spectrum at all energies can be well approximated by the 
simple point source solution. 



5 APPLICATION TO THE NEARBY SNRS 

In this section, we shall apply our study to the nearby known SNRs 
listed in Table 1 with distances < 1 kpc from the Earth. It should be 
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Table 1. Parameters of known SNRs with distances < 1 kpc from tiie Earth. 
References: (1) Blair et al. 2005 (2) Tatematsu et al. 1990 (3) Leahy & 
Aschenbach 1996 (4) Leahy & Tian 2007 (5) Braun et al. 1989 (6) Caraveo 
et al. 2001 (7) MiceU et al. 2008 (8) Jian et al. 2005 (9) Strom 1994 (10) 
Plucinsky et al. 1996 (11) Caraveo et al. 1996 (12) Egger & Aschenbach 
1995 (13) Katsuda et al. 2008 (where for the age we take the mean of 2700 
yr and 4300 yr reported in the Uterature). 



SNR 


Distance (kpc) 


Age (yr) 


References 


Cygnus Loop 


0.540 


1.0 X lO-* 


1 


HB21 


0.800 


1.9 X 10* 


2,3 


HB9 


0.800 


6.6 X 103 


4 


S147 


0.800 


4.6 X 10^ 


5 


Vela 


0.294 


1.12 X 10-* 


6,7 


G299.2-2.9 


0.500 


5.0 X 10^ 


8 


SN185 


0.950 


1.8 X 10^ 


9 


Monogem 


0.300 


8.6 X 10"* 


10 


Geminga 


0.157 


3.4 X 10"' 


11 


Loopl 


0.170 


2.0 X lO"' 


12 


Gl 14.3+0.3 


0.700 


4.1 X 10"' 


8 


Vela Junior 


0.750 


3.5 X 103 


13 



mentioned that some of the age and the distance parameters listed 
in Table 1 carry large uncertainties. For instance, the distance to 
the Geminga was measured to be 157 pc using Hubble Space Tele- 
scope (HST) observations (Caraveo et al. 1996) but recently, again 
using HST measurements, Faherty et al. 2007 reported the distance 
of Geminga to be 250 pc. For the Cygnus Loop, Minkowski 1958 
reported a distance of 770 pc whereas measurements based on HST 
observations claimed a distance of 440 pc (Blair et al. 1999). Re- 
cent measurements further claimed the distance to be 540 pc (Blair 
et al. 2005). For HB21, Tatematsu et al. 1990 measured a distance 
of 800 pc and Leahy & Aschenbach 1996 estimated an age of 
1.9 X 10* yr while later, Byun et al. 2006 suggested a distance of 
1.7 kpc and Lazendic & Slane 2006 estimated an age of 5.6 x 10^ 
yr. Leahy &: Aschenbach 1995 estimated the distance and age of 
HB9 as 1 kpc and 7.7 x 10'' yr respectively, and Leahy & Tian 
2007 suggested a distance of 800 pc with sedov age of 6.6 x 10^ 
yr and age of (4 — 7) x 10^ yr based on evaporation cloud model. 
The lack of precise informations on these parameters can affect our 
results because of the strong dependence of the CR spectrum on 
these parameters. 

For our study, we will assume that the proton source index 
r = 2.13 so that for 5 = 0.6, we get F -h 5 = 2.73 the observed 
proton spectral index at the Earth (Haino et al. 2004). It should be 
noted that the value of the source index can depend on the choice of 
the propagation model and different propagation models may take 
different values. For instance, models based on diffusive reaccel- 
eration in the Galaxy favors a diffusion index of 5 ~ 0.3 which 
corresponds to a source index of F ~ 2.4 (Trotta et al. 2011). 
This is steeper than the value adopted in our present work which 
is based on a purely diffusive model of CR propagation. For the 
CR electrons, to get the source index, we first determine the back- 
ground spectrum. This is done by fitting the observed data between 
(10 - 200) GeV provided by the FERMI and the PAMELA ex- 
periments (Ackermann et al. 2010b, Adriani et al. 2011). We as- 
sume that this is the energy region where the contamination due 
to the local sources as well as the effect of the solar modulation 
are minimum. From the fit, the background spectral index is found 
to be 3.10 ± 0.01. Under diffusive propagation model, CR elec- 
trons produced by a uniform and stationary source distribution, and 
subject to radiative losses during their propagation in the Galaxy 
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Figure 9. CR proton spectra from nearby SNRs listed in Table 1 for the 
three different source models: point source (solid line), surface source 
(dashed line) and solid source (dotted line). The contributions from the in- 
dividual SNRs are labelled by their names and their total contributions as 
"Total local SNRs". The thick sohd Une represents the fitted total observed 
spectrum, 1.37 x (B/GeV)-^ ''^ cm-^ s"! sr"! GeV"! taken from 
Haino et al. 2004. Our calculation assumes an energy independent burst- 
like injection of particles at ta = and that each SNR produces 10®'' ergs 
of CR protons. Other model parameters: V = 2.13, Do = 2.9 x 10^* cm 
2 s-i, 5 = 0.6 and Eo = 3 GeV. 
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Figure 10. Same as in Fig. 9 but for an energy dependent escape model. 
We assume a = 2.0, Emax = 10® GeV and t^ed = 500 yr. The dashed 
line represents the "Total local SNRs" we obtained in the case of the point 
source approximation shown in Fig. 9. 

results into an equilibrium spectrum given by where 
/3 = (1 - S)/2 (see e.g. Thoudam & Horandel 2011). Using the 
value of the background index obtained from the fit, we get the 
electron source index as F = 2.3. This is the value we will adopt 
for the rest of our calculations for the electrons. Furthermore, in 
the following we will assume that 10% of the supernova explosion 
energy of 10®^ ergs converts into CR protons and 0.1% into the 
electrons. All these parameters are assumed to be the same for all 
the SNRs. 



5.1 CR protons 

First, we compare the results of the point source approximation 
with those of the spherical solid and the surface source models 
which are all based on the energy independent escape model. These 
are shown in Fig. 9 where different lines represent different source 
model: thin solid (point), dashed (surface) and dotted (solid). The 
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thick solid line represents the fitted observed proton spectrum given 
in Haino et al. 2004. The calculation assumes a burst-like injec- 
tion of CRs at time to = 0. The contributions from the dominant 
SNRs in different energy intervals are indicated by their names. Al- 
though some of the individual SNRs like the Vela and the G299.2- 
2.9 show different spectra at low energies under the different source 
models, the differences are not significant in the total combined 
spectrum from the nearby SNRs. It is because at low energies be- 
low ~ 60 GeV, the dominant contributions are from the Mono- 
gem and the Loopl whose spectra does not show any differences 
between the models because of their old ages. Although different 
SNRs dominate at different energy intervals, as a whole their total 
spectrum looks smooth except for a slight dent somewhere between 
~ (10 - 100) GeV. 

We have also checked the results for other values of to ^ 10^ 
yr. Except for the signatures of the absence of young sources as 
to takes larger values, we have found that their total spectra does 
not show any significant differences between the different source 
models at all values of to- 

When we apply the energy dependent escape model, we find 
that the results are significantly different from those of the energy 
independent models. This is shown in Fig. 10. The dashed line rep- 
resents the total spectrum we obtain for the point source approxi- 
mation shown in Fig.9. We can see that the total spectrum in the 
energy dependent case show irregular structures which are due to 
the low energy spectral breaks of the individual SNRs. However, 
such features can possibly remain embedded in the dominant CR 
background and may not be distinctly visible in the observed spec- 
trum. 

5.2 Electrons 

The electron spectra for the point, solid and the surface source mod- 
els are shown in Fig. 11 . In the figure the data are from the FERMI 
(Ackermann et al. 2010b), PAMELA (Adriani et al. 2011) and the 
HESS (Aharonian et al. 2008c, 2009) experiments. As in the case 
of protons, the total electron spectra also does not show any dif- 
ferences between the different models. However, unlike in the case 
of the protons, the total electron spectra show some irregular fea- 
tures near the highest energies which are due to the effects of sharp 
cut-offs in the individual spectra due to radiative energy losses. 
For instance, the strong peak at iS ~ 10^ yr is due to the effect 
of SN185. From the figure, we can notice that at energies greater 
than few TeVs our results which are based on a pure power-law 
source spectrum significantly over predicts the data. Taking larger 
values of to can suppress the contributions of Vela, G299.2-2.9 and 
SN185 which are the dominant contributors at high energies. For 
to = 2 X 10* yr, their contributions will be completely removed. 
This points towards the importance of source modeling in order to 
understand the contribution of local sources in the high energy elec- 
tron spectrum. One common way to handle this problem is to as- 
simie an exponential cut-off exp{—E/Ec) in the source spectrum 
at a few TeVs (see e.g. Delahaye et al. 2010). Another possibility 
is that the high energy electrons might have suffered significant en- 
ergy losses within the SNR itself before they are released into the 
ISM (Thoudam & Horandel 2011). Therefore, electrons at higher 
energies might be released with a spectrum steeper than the lower 
energy ones. For the present study, we adopt the much simpler ex- 
ponential cut-off and in Fig. 12, we show the results obtained for 
Ec = 2 TeV. We can see that the shape of the total local spectrum 
in the TeV region is now determined mostly by the exponential cut- 
off and the irregular structures present near the highest energies in 
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Figure 11. Electron spectrum from the nearby SNRs listed in Table 1 for 
the three different source models: point source (solid line), surface source 
(dashed line) and solid source (dotted Une). We assume a burst-like injection 
of particles and a pure power-law source spectrum of index T = 2.3 with 
each SNR producing 10*^ ergs of CR electrons. AH other model parameters 
remain the same as in Fig. 9. The data are taken from the FERMI, PAMELA 
and HESS experiments. 
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Figure 12. Same as in Fig. 1 1 but for a source spectrum with an exponential 
cut-off at Ec = 2 TeV. The thick dashed line represents the background 
spectrum (see text for details) and the thick solid line represents the total 
background plus nearby SNRs in the point source approximation. The data 
are the same as in Fig. 11. 

Fig. 1 1 no longer exist. In Fig. 12, the thick dashed line represents 
the background spectrum (which we obtain as mentioned before) 
with an index 3.1 and an exponential cut-off at 2 TeV. The thick 
solid line represents the total background plus the nearby SNRs 
contribution obtained in the point source approximation. Detailed 
calculations of the background spectrum taking into account the 
various source models discussed here will be presented elsewhere. 

For the energy dependent escape model, the results are shown 
in Fig. 13 for a pure power-law source spectrvun. In the figure, we 
also show for comparison the total local spectra obtained in the 
case of the point source approximation (dashed line in Fig. 11). The 
total spectrum show several irregular features and spikes. These 
features are stronger than the ones present in the proton spectrum 
which is due to the presence of additional breaks in the individual 
electron spectra at high energies. The position of these spikes not 
only depends on the age and distance of the individual SNRs but 
also on the assimie energy dependent escape model (i.e., on the 
parameters a, tsed and Emax)- 

In Fig. 14, we show the electron spectra obtained using the 



12 Satyendra Thoudam 



Point source 

PAMELA 
HESS 2008 
HESS 2009 




10^ 10^ 10^ 10* 10^ 



E (GeV) 

Figure 13. Electron spectrum from nearby SNRs for the energy dependent 
injection model and a pure power-law source spectrum. The dashed line 
represents the "Total local SNRs" obtained under the point source approx- 
imation shown in Fig. 11. We assume a = 2.0, Emax = 10® GeV and 
*sed = 500 yr. The data and all other model parameters remain the same 
as in Fig. 11. 
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Figure 14. Same as in Fig. 13 but for a source spectrum with an exponential 
cut-off at Ec = 2 TeV. The thick dashed Une represents the background 
spectrum and the thin dashed line represents the "Total local SNRs" for the 
point source approximation shown in Fig. 12. The data are the same as given 
in Fig. 11. 

same model as in Fig. 13 but with an exponential cut-off in the 
source spectrum at Ec = 2 TeV. In the figure, the thin and the 
thick dashed lines represent the local SNRs contribution in the case 
of the point source approximation and the background spectrum 
respectively as shown in Fig. 12. The thick solid line denotes the 
total background plus local spectrum for this case. We can notice 
that even after imposing the cut-off some prominent features still 
remain at few TeVs in the overall total spectrum imlike in the point 
source approximation where the cut-off almost smoothens the total 
spectrum. 



6 OVERALL RESULTS AND DISCUSSIONS 

We have shown that the commonly adopted point source approx- 
imation does not always remain a good approximation as far as 
nearby CR sources are concerned. For a typical source distance of 
Ts ~ (100 — 300) pc, we have shown that at low energies, the point 
source results for young sources (t3'pically t < 10* yr) differ sig- 
nificantly from those calculated using a finite source size. At high 



energies the point source still remains a good approximation. Un- 
der the energy independent particle escape model, we found that 
the effects of the finite source size are similar in both the types of 
particle injection model considered in our study: the burst-like and 
the continuous injection. For very old nearby sources it > 10 ' yr), 
we have found that the results are independent of both the source 
size and the particle injection model and hence, the burst-like point 
source model represents a good approximation at all energies. We 
have also shown in Fig. 3 that for a given value of the CR diffu- 
sion coefficient, there is a certain parameter space in (vs, t) under 
which the point source approximation remains valid for CRs of our 
interest, i.e., with energies S > (3 — 10) GeV. When applied to 
the nearby known SNRs within 1 kpc, interestingly we have found 
that their total spectrum almost remain the same in the three differ- 
ent source models although some of the individual SNRs like Vela 
show differences between the models (Figs. 9 & 11). We found that 
it is because at low energies where the point source approximation 
is most likely to break down, the local spectrum is dominated by the 
Monogem and the Loopl. These SNRs are quite old with Monogem 
age 8.6 x 10'* yr and Loopl ~ 2.0 x 10* yr due to which their 
CR spectra at the Earth are independent of their sizes and are well 
represented by the point source solutions. 

We have also studied an energy dependent escape scenario 
where CRs of different energies are assumed to escape at differ- 
ent times during the SNR evolution. We assumed that the escape 
time follows, tesc oc E~^^°' with a chosen to be equal to 2.0. 
Under this model, we assumed that the highest energy particles es- 
cape the remnant at the start of the Sedov phase followed by the 
lower energy ones at later times. For Emax = 1 PeV, tsed = 500 
yr and the maximum CR confinement time of 10* yr adopted for 
our study, we found tesc = (500 — 10*) yr and the escape ra- 
dius Resc = (5 - 100) pc for energies E = (IPeV - 25GeV). 
For young sources, the spectrum obtained under this model show 
breaks at lower energies which are due to the longer confinement 
times at those energies. At high energies, the results are very similar 
to those of the point source approximation. This is not just because 
of the small values of Resc at high energies but also due to their 
large values of D{E) at these energies. In fact, we have shown in 
section 2 that even for a large escape radius of 100 pc, the point 
source still represents a good approximation at high energies (see 
e.g., Fig. 1 surface source). Therefore, it should be understood that 
it is not the small Resc, but actually the large D{E) which is re- 
sponsible for the point source validity at high energies under the en- 
ergy dependent escape model. When applied to the nearby known 
SNRs, we have found that the results obtained under this model are 
significantly different from those obtained under the point source 
approximation. The total local spectrum show more irregular struc- 
tures as compared to the point source results. Also, we have noticed 
that there is a big dip between around (10^ — 3 x 10^) GeV which 
is mainly due to the low energy cut-off in the Vela spectrum (Figs. 
10 & 13). These results seem to suggest that if SNRs are the main 
sources of CRs in our Galaxy, then the widely adopted point source 
approximation with an energy independent escape scenario appears 
flawed for CR studies from the nearby SNRs. 

For the protons, the irregular spectral features that we have 
found in the energy dependent escape model may be suppressed by 
the dominant background produced by distant sources and hence, 
may not show up distinctly in the total observed spectrum. But, 
for the electrons they can possibly show up to detectable levels es- 
pecially at TeV energies where the background level is expected 
to be significantly less. Recently, Kawanaka et al. 2011 proposed 
that such spectral features can be used to estimate the CR con- 



On the point source approximation of nearby cosmic-ray sources 13 









8 






/ fj 















E (GeV) 



E.= 20TeV - 
E^= 5 TeV - 
Er= 2 TeV .. 




E (QeV) 

Figure 15. Top: Total electron spectrum from nearby SNRs assuming a 
pure power-law source spectrum [Ec = oo) for different values of a: 1.8 
(solid line), 2.0 (dashed line), 2.5 (dotted line). Bottom: For a = 1.8 with 
exponential cut-off in the source spectrum at Ec = 10 TeV (solid Une), 
2 TeV (dashed Une) and 1 TeV (dotted Une). The data and other model 
parameters remain the same as in Fig. 11. 



finement time inside SNRs. Their study assumed a single nearby 
source having characteristics similar to that of the Vela remnant. 
It should be noted that the position and the strength of such fea- 
tures strongly depend on the CR escape model especially on the 
a parameter and also on Ec (if there is an exponential cut-off in 
the source spectrum). For instance, assuming Ec > 2 TeV would 
produce stronger features and vice versa compared to our results 
shown in Fig. 14. Similarly, assuming a < 2.0 would produce 
stronger peaks at comparatively higher energies as low energy CRs 
would be confined for relatively longer times as indicated by Figs. 
6 & 7, and taking q > 2.0 would smoothen the peaks as low en- 
ergy CRs would also start escaping at early times. For a » 2.0, 
the energy dependent results will tend towards the point source re- 
sults obtained for to = tsed- These can be understood from Fig. 
15 where in the top panel we have shown the electron spectrum for 
the case of pure power-law source spectrum (which corresponds to 
Ec = oo) for different values of a: 1.8 (solid line), 2.0 (dashed 
line), 2.5 (dotted line). We can clearly see the left shifting of the 
peak between (1 — 10) TeV as a increases from 1.8 to 2.5. In the 
bottom panel, we have shown the spectra calculated for the case of 
a = 1.8 for source spectra with exponential cut-offs at Ec — 10 
TeV (solid line), 2 TeV (dashed line) and 1 TeV (dotted line). Here 
again, we can notice that the peak at around 5 TeV grows stronger 
as Ec takes larger values. 

If we look into the HESS electron data, there is an indication 
of an abrupt rise at the highest measured energy. If future better sen- 
sitive experiments like the CTA and the CALET provide good qual- 



ity data at these energies, that would indeed provide useful infor- 
inations to understand CR escapes from some of our nearby SNRs. 
However, the large uncertainties involved in the age and the dis- 
tance estimates of some of these SNRs may be an issue because of 
the strong dependence of the CR spectrum on these parameters. Re- 
garding this, measurements of electron anisotropy both amplitude 
as well as its direction at these energies might also be important in 
order to identify the dominant source. 

Recently, Di Bernardo et al. 2011 studied the contributions 
of the nearby pulsars and the SNRs to the high energy electron 
spectrum. One of their conclusions is that a strong contribution 
from the nearby SNRs is not supported by the recent upper lim- 
its on the electron anisotropies provided by the FERMI Large Area 
Telescope observations (Ackermann et al. 2010a). But, it should be 
noted that their calculations assumed the sources to be burst-like 
point sources emitting CR particles independent of energy. In Fig. 
14, we show that between ~ (lOOGeV — 2TeV), the contribution 
from the nearby SNRs is significantly larger in the point source ap- 
proximation than in the energy dependent model. We believe that 
a more realistic treatment of particle escape model from the SNRs 
may change their conclusion. Other class of sources which might 
also produce significant contributions to the high energy leptonic 
(electron plus positron) spectrum are pulsars and dark matter. Mod- 
els based on these sources are motivated mostly by the detection of 
the rise in the positron fraction above ~ 10 GeV by the RAMELA 
experiment (Adrian! et al. 2009). If we assume that positrons are 
produced only during the interaction of the priinary CRs with the 
interstellar gas, the positron fraction is expected to decrease with 
energy which is in contrast to the observations. A possible solu- 
tion to this problem among others may be the presence of one or 
more nearby positron sources like pulsars or dark matter (see e.g. 
Grasso et al. 2009 and references therein). Future measurements of 
electron anisotropies with better sensitivities and also the absolute 
positron spectrum at high energies can provide better understand- 
ing of the nature and the type of the dominant source(s). 

Moreover, a good understanding of the background contribu- 
tion would also be crucial. In an earlier paper, we had presented 
calculations of the averaged background based on a simple energy 
independent model of CR confinement within the SNRs (Thoudam 
& Horandel 201 1). In future, we will present background estimates 
for both the protons and the electrons taking into account the energy 
dependent confinement/escape of particles. The calculation will in- 
clude the various energy loss and the interaction processes taking 
place during the time particles are confined within the sources. In 
addition, we will also present the possible effects on other observed 
CR properties like the Galactic diffuse 7-ray emission, s/p ratios 
and the anisotropies. 
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